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Actin is a globular protein which forms long polar filaments in the eukaryotic cytoskeleton. Actin 
networks play a key role in cell mechanics and cell motility. They have also been implicated in 
information transmission and processing, memory and learning in neuronal cells. The acting fila¬ 
ments have been shown to support propagation of voltage pulses. Here we apply a coupled nonlinear 
transmission line model of actin filaments to study interactions between voltage pulses. By assigning 
a logical Truth to the presence of a voltage pulses in a given location of the actin filament, and 
False to the pulse’s absence we represent digital information transmission along these filaments. 
When two pulses, representing Boolean values of input variables, interact, then they can facilitate 
or inhibit further propagation of each other. We explore this phenomenon to construct Boolean 
logical gates and a one-bit half-adder with interacting voltage pulses. We discuss implications of 
these findings on cellular process and technological applications. 


Actin is a globular structural protein; one of the most 
highly conserved eukaryotic protein, ranging from uni¬ 
cellular organisms to plants and animals. Actin plays a 
key role in cell motility forming actin filaments (microfil¬ 
aments), which in turn generate parallel bundles that the 
cell utilizes for contractility, in cell division for the cre¬ 
ation of the cleavage furrow and in cell motility. Actin 
has also been correlated with nervous system activity and 
learning: actin cytoskeleton affects synaptic properties 
leading to learning, actin and its regulatory proteins are 
involved in various stages of memory [lH5]. The actin fil¬ 
ament networks are key players in modulation of synaptic 
terminals due to their interactions with ion channels [6j 
and filtration of noise in synapses cub]- By modulating 
dendritic ion channel activity actin filaments govern neu¬ 
ral information processing and facilitate computational 
abilities of dendritic trees via facilitation of ionic con¬ 
densation and ion cloud propagation [9]. 

Previously we proposed a model of actin filaments as 
two chains of one-dimensional binary-state semi-totalistic 
automaton arrays to signalling events, and discovered lo¬ 
cal activity rules that supports travelling or stationary 
localizations m- This finite state machine model has 
been further extended to a quantum cellular automata 
(QCA) model in [JTj. We have shown that quantum actin 
automata can perform basic operations of Boolean logic, 
and implemented a binary adder m, and three valued 
logic operations [12] . These models were implemented in 
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general algorithmic terms, without describing any spe¬ 
cific physical mechanisms that could be used to imple¬ 
ment cell excitations and interactions. They also did not 
employ interactions between propagating localisations in 
a spirit of collision-based computing fl3j . 

We aim to rectify these omissions and present a model 
of actin in terms of RLC (resistance, capacitance, induc¬ 
tance) non-linear electrical transmission wires that can 
implement logical gates via interacting voltage impulses. 
These models take inspiration from the idea of a cellu¬ 
lar automaton, but do not fit exactly in its definition, as 
cells do not possess discrete states, there are no intrinsic 
time steps, and no rules are defined for state transitions. 
However, the states can be digitised by defining a suit¬ 
able threshold, a time step can be defined by convenience 
and transitions are governed by Kirchhoff’s circuit laws 
so that a deterministic evolution takes place even in ab¬ 
sence of explicit transition rules. 

Our starting point is the usual model of electrical wires 
as sequences of circuits composed by resistors, capaci¬ 
tors and inductors. Our model is based on Tuszynski et 
al. HU model of actin monomers in terms of electrical 
components. This latter model was developed in order 
to explain experimental observations of ionic conductiv¬ 
ity along actin filaments m • This model exhibits, in a 
continuous limit, solitons (standing waves and travelling 
impulses). 

The aim of the present work is to use the same equa¬ 
tions as in. [14j but without invoking the continuous 
limit approximation and to study the behaviour of sev¬ 
eral types of solutions for tens of monomers. We will 
show how an excitation moves along the actin filament 
and how it collides with another excitation coming from 
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elsewhere; we will show that these collisions can be used 
to implement logical operations. This could provide a 
physical description of signal propagation and processing 
in the cellular milieu and also possibly in hybrid bio-nano 
technological applications. 


I. ACTIN FILAMENTS AS NONLINEAR RLC 
TRANSMISSION LINES. 



FIG. 1. A circuit diagram for the n-th unit of an actin filament 
(from [H]). 


In this section we recall the basic model that will be 
used throughout the paper; more details can be found in 
the original paper M- 

Referring to Fig. [l] where an actin monomer unit in a 
filament is delimited by the dotted lines, we assume that 
capacitors are nonlinear (see discussion in m, El): 


Qn = Co(V n - bV%) (1) 

The main equation is: 

LC 0 ^(V n -bVZ) = 

= V n+1 + V-i - 2V - i?iC 0 |(V - bV%) 

- R 2 C 0 {2f t (V n - &V 2 ) - i(V n+ 1 - &V 2 +i) 


where e is the electrical charge, cq the permittivity of 
the vacuum, e the dielectric constant of the solution 
where the actin molecule is placed, estimated similar to 
e water ~ 80, ks Boltzmann’s constant and T the absolute 
temperature, that is taken « 293 K. The Bjerrum length 
is the distance beyond which thermal fluctuations are 
stronger that electrostatic interactions between charges. 
The parameter used in m is A b = 7.1 x 10 10 m. 

The formula for the capacitance is as follows: 


Co 


27reeo / 


ln( 


T_actin~ l ~\b n 

R actin J 


( 4 ) 


where l is the length of a monomer, taken ~ 5.4nra; 
and r act in ~ 2.5 nm. With the parameters given, the 
capacitance per monomer is estimated Co « 96x 10 ~ 6 pF. 
For the inductance we have: 

j _ /JjN 7t(t actin T A_g) . . 


where p is the magnetic permeability and N is the num¬ 
ber of turns of the coil, that is the number of windings of 
the distribution of ions around the filament. It is approx¬ 
imated by counting how many ions can be lined up along 
the length of a monomer as N = l/r and it is supposed 
that the size of a typical ion is ~ 3.6 x 10 -10 m. With 
these assumptions we have for a monomer L « 1.7 pH. 

For the resistance we have 


R = 


p /?7/( {v ac tin T / Tactin) 

2irl 


( 6 ) 


where resistivity p is approximately given by: 


A o + c K + + ^ o a + C Na + 

where c K + and c Na + are the concentrations of sodium 
and potassium ions, considered respectively 0.15 M 
and 0.02 M; A^ + « 7.4(Om) _1 M _1 and Aff a+ ~ 
5.0depends only on the type of salts. Ac¬ 
cordingly, Ri « 6.11 MQ; R 2 is taken = Ri/7 « 0.9 MQ 
following [14]. 

We note that resistance, inductance and capacity val¬ 
ues could be tuned by changing the temperature and the 
ion concentrations. Further sections deal with solutions 
of © and their interactions. 


- |(V-i-&V 2 -i)l (2) 

The above is the basic equation, that will be used to 
compute the cells’ behaviour. 

In p3] formulas are derived to obtain values for the 
relevant parameters. The Bjerrum length is defined as 
follows: 

47reeo ksT ® 


II. EVOLUTION OF PULSES 

We recap the parameters that will be used in almost 
all computations: 

L = 1.7p H, Co = 96 • 10- 6 pF, ( . 

R 1 = 6.11 • 10 6 U, R 2 = 0.9 • 10 6 U w 

With these parameters we measure time in nanoseconds. 
Let us consider the case when an initial voltage is set 
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FIG. 2. An input pulse travelling along act in filament: (a) pa¬ 
rameters (b) the same parameters but Ri = 61.110 5 O. 
The pictures represent the cell evolutions, each cell voltage 
is represented by a diagram with time on the horizontal axis 
and voltage V on the vertical one. The diagrams are tiled 
vertically with cell 1 at the bottom and cell 20 at the top, 
with a unitary displacement. 

at one end of the filament for some time as, for exam¬ 
ple, could take place when an act in filament is in close 
proximity of an ion channel. We use for voltage the form: 

! to _g t~\~to 

V ° ~ 2 2(e t ~ t ° + e~ t+t °) 9 

In all the following computations we consider to = 3 ns. 

We find that pulses actually travel, but after a few 
nanoseconds they are damped (Fig. IT)- From the dia¬ 
gram we see that the 12th unit is reached in about 4 ns, 
so that the speed is approximately 5.4 nm x 12/4 ns, that 
is about 16 m/sec. That is consistent with the estimates 
in [14] . If we decrease the resistance R\ by an order of 
magnitude, the speed doubles. Then we can obtain a 
much longer lasting pulse (Fig. IT)- 

Figure [3^i shows the case where two identical impulses 
are applied at the ends of the act in filament. In Fig. ^jp 
we see interaction of two pulses of opposite sign applied 
to ends of the filament. 


FIG. 3. Two input pulses travelling along actin monomers 
from the opposite ends, with parameters (|8|, but R± = 
61.110 5 U. (a) Identical pulses, (b) Opposite voltage. 

In the next two sections we consider two types of in¬ 
puts: forced and unforced. An unforced input is an im¬ 
pulse applied just once at input sites of actin-based gate. 
A forced input is an impulse applied continuously during 
evolution of the gate. 

III. GATES WITH UNFORCED PULSES 

Let us consider interactions and possible gates imple¬ 
mented using only initial conditions of cells, without any 
forced pulses applied. They are based on the assumption 
that we can initialize any cells with a pre-determined 
electrical potential value and to reliably measure output 
potentials at given units of actin filament. 

Using parameters ^ and b = 0.1, and initial condition 
V = Vo 7^ 0, we look at cells having voltages greater than 
a fixed percentage of initial potential Vo- We find that 
some logical operations are possible to implement, if we 
consider signals above « 0.1 Vo as output, to read outputs 
of the logical gates. 

Figure [4] shows the AND operation obtained with two 
input cells at positions 8 and 15. The output is in cells 11 
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FIG. 4. Unforced pulses, AND operation. Inputs are (a) 00, 
(b) 10, (c) 01 and (d) 11 in cells 8 and 15; output in cells 11 
and 12. Cells excited above 0.1 Vo are black; time evolves along 
the horizontal axis and cells are aligned along the vertical axis. 


FIG. 5. Unforced pulses, OR operation. Inputs are (a) 00, 
(b) 10, (c) 01 and (d) 11 in cells 9 and 11; output in cell 10. 
Cells excited above 0.1 Vo are black; time evolves along the 
horizontal axis and cells are aligned along the vertical axis. 


and 12, that are excited at the end of the evolution if both 
the input cells were excited initially. Figure [5] shows the 
OR operation obtained with two input cells at positions 
9 and 11. The output is in cell 10 that is excited at the 
end of the evolution if one of the input cells or both of 
them were excited initially. 

To obtain a not gate, we initialize an auxiliary cell at 
position 9 at —Vo and place the input cell at position 11; 
the output is found in cell 10 as shown in Fig. [6] 

To build an XOR gate we use the formula 

a ® b = (-i(a A b )) A (a V b) (10) 

as shown in Fig. [7] Each run is composed by four gates, 
and it is supposed that some mechanism reads the output 
of each gate, amplifies the signal and passes the signal to 
the next gate. In the figure the arrays represent signals 
travelling between gates. The first from the left is an AND 
gate with input cells 8, 15 and output cell 12; its output 
goes to the next gate: the NOT gate having input in cell 
11 and output in cell 10. The original inputs, coming out 
at cells 8 and 15 of the first gate, are sent to input cells 
9 and 11 of the OR gate in the third section. Finally, the 
output of the OR gate in cell 10 and the output of the 
not gate are sent to input cells 8 and 15 of the final AND 
gate. The final output of the XOR operation is found in 
cells 11 and 12 at the right end. 



0 5 10 15 20 25 30 0 5 10 15 20 25 30 


(a) (b) 

FIG. 6. Unforced pulses, not operation. Input is (a) 0 and 
(b) 1 in cell 11; output in cell 10. Cells excited above 0.1 Vo 
are black; time evolves along the horizontal axis and cells are 
aligned along the vertical axis. 


IV. GATES WITH FORCED PULSES 

Let us we consider interactions and possible gates im¬ 
plemented by applying input pulses continuously at input 
sites. Namely, we send a sinusoidal pulse to input during 
the whole evolution of the system, with the intensity of 
the input pulse remaining constant. 

We use slightly different values for resistances: R\ = 
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FIG. 7. Unforced pulses, XOR operation. Inputs are (a) 00, 
(b) 10; (c) 01 and (d) 11 in cells 8 and 15 of the leftmost 
gate; output in cells 11 and 12 of the rightmost. Cells excited 
above 0.1 Vo are black; time evolves along the horizontal axis 
and cells are aligned along the vertical axis. 



FIG. 8. Forced pulses, and operation. Inputs are (a) 00, 
(b) 10, (c) 01 and (d) 11 in cells 8 and 14; output in cell 11. 


9.23 10 6 £2, i7 2 = 1.3210 6 £2, that can be obtained chang¬ 
ing concentrations of K + and 7Va + ; excitations are sinu¬ 
soidal with a 1 ns period. 

Figure [8] shows the AND operation obtained with two 
input cells at positions 8 and 14. The output is found in 
cells 11 that is excited at the end of the evolution if both 
the input cells are excited. 

Figure [9] shows the XOR operation obtained with two 
input cells at positions 10 and 14. The output is found in 
cells 12 that is excited at the end of the evolution if only 
one the input cells is. The input pulses are sinusoidal 
with opposite phases. The other logical gates can be 
implemented in a similar way. 

We also show in Fig.[l0]a one-bit half-adder: it is com¬ 
posed by an AND gate followed by an XOR. The AND gate 
has input cells 8 and 14 and outputs the most significant 
bit (the carrier) in cell 11; the XOR has input cells 10 and 
14 and outputs the less significant bit in cell 12. It is 
supposed that input in cell 14 is common to both gates, 
while input in cell 8 of the AND gate travels to cell 10 of 
the XOR changing its phase. 

One could argue that it could be difficult in practice 
to send an input pulse to exactly one monomer (or to 
initialize it in an excited status), which requires a local¬ 
ization precision of under 5 nm. We therefore compute 
some evolutions supposing that the pulse hits a number 
N of monomers. Under these conditions, the values of re¬ 
sistance, inductance and capacitance must be computed 



(a) (b) 




(c) (d) 

FIG. 9. Forced pulses, XOR operation. Inputs are (a) 00, 
(b) 10, (c) 01 and (d) 11 in cells 10 and 14; output in cell 12. 
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with appropriate addition rules. For resistance we have: 

N N 

R2tot = (^2 ^2,i ) \ Rltat = (11) 


Total inductance and capacitance are obtained sum¬ 
ming those of individual monomers. Figure [TT] shows the 
XOR gate with input sent to groups of 10 monomers in 
positions 8-17 and and 23-32. As in Fig. [9j the input 
pulses are sinusoidal with opposite phases. Output is in 
cell 20. In our approximation, we do not consider the de¬ 
tailed evolution inside the groups of monomers receiving 
the input, they are treated as a single cell, with R, L, C 
parameters computed summing on 10 monomers. 

We note that the results are robust: pulses fade slowly 
and a threshold of 0.3 of the input signal can be used, 
instead of 0.1 as in the previous computations. 


V. CONCLUSION 

In numerical experiments presented in this paper we 
have shown that it is possible to build basic logical gates 
from actin filaments considered as RLC circuits, provided 
that input pulses can be sent to specific locations and 
initial values of voltage can be set in actin monomers. 
We demonstrated that several logical gates can be imple¬ 
mented via interactions between input pulses. For exam¬ 
ple, to realise an AND gate we send the same phase pulses 
to the input cells. The pulses sum with each other. The 
output threshold is set in such manner that if there is 


.1,R1=9.23*10 A 6,R2=1.32*1G A 6,L=1.7*10 A -12,CO=96*10 A 


,1,R1=9.23*10 A 6,R2=1.32*10 A 6,L=1.7*10 A -12,CO=96*10 A 
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FIG. 11. Forced pulses, XOR operation. Inputs are (a) 00, 
(b) 10, (c) 01 and (d) 11 in cells 8-17 and 23-32; output in 
cell 20. 


only one active input the output potential is lower than 
the threshold. If there are two summed input potentials 
the threshold is exceeded. To realise an XOR gate we send 
pulses with opposite phases so when they meet they are 
subtracted from each other. 

To try a real implementation it is necessary to specify 
how output can be measured with sufficient sensitivity 
(either 0.1 or 0.3 of the initial voltage) and how output 
pulses of a gate can be reliably sent as input to another. 
These points will be considered in future works, we can 
however anticipate some considerations about cascading 
gates. 

Actin filaments are usually found in networks, form¬ 
ing cross links that can have a parallel or orthogonal 
shape. These links can be made by several kinds of bind¬ 
ing proteins that originate different network topologies. 
We consider for instance filamin A (see e.g. [18]), that 
builds orthogonal networks: where two actin filaments 
cross link, we suppose that electric pulses generated as 
output of a gate in a filament can reach the other using 
the filamin molecule as a connecting wire. Branching of 
actin filaments networks has been observed in cells and 
it requires the presence of ARP2/3 protein as a branch 
point from which two actin filaments propagate at acute 
angles m- 

Considering for the sake of simplicity the same electri¬ 
cal mechanisms as for actin, we could compute Ri, 

Co, and L for filamin, taking into account just its dimen- 
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sions. 

In this way, in principle, any number of actin gates 
could be connected, however there is a problem: pulses 
intensity will fade away rather quickly, unless some am¬ 
plifying method is found. 

Some studies (e.g. [20] ) have found additional roles of 
filamin A in signal transduction, and others like [2T] have 
found that actin filaments can regulate some voltage- 
gates ion channels. We can therefore think of a system 
where the output voltage of a gate opens a ion channel 
that in turn activates the next gate and so on. Of course 
this subject must be treated in more depth to draw any 
practical conclusions, as the spatial distribution of the 
ion channels within the actin network, their dimensions, 
the necessary electrochemical gradients, the type of ions 
that can be used, their diffusion rates etc. must be care¬ 
fully evaluated to check the feasibility of a logical device. 

Tubulin is another important structural protein that 
polymerises into cylindrical microtubule structures that 


have been found to be electrical amplifiers [22]. Actin fil¬ 
aments are connected to tubuline microtubules via micro¬ 
tubule associate proteins (MAP) [23, 124] , Thus, if elec¬ 
trical properties of the MAP are similar to that of actin 
units we can speculate that cascading of Boolean gates 
into arithmetic circuits can be done via actin to MAP 
to tubuline to MAP to actin links. It has been recently 
proposed that microtubule phosphorylation by calmod¬ 
ulin kinase II enzyme complex can be linked to memory 
formation at a cellular level [25] . In a manner similar to 
that presented in this paper for actin filaments, signals 
propagating along microtubules have been also shown to 
be capable of interacting via logical gates [25]. This can 
lead to an integration of signal storage and processing ca¬ 
pabilities of the cytoskeleton, which is composed of both 
actin filaments and microtubules that interact with each 
other. 

In case of in vitro desing of actin-based computing cir¬ 
cuits a tailored mono-molecular amplifiers [26 can be 
inserted between actin strands in cascaded gates. 


[1] J. Tuszynski, J. Brown, P. Hawrylak, and P. Marcer, 
Philosophical Transactions Royal Society of London. Se¬ 
ries A: Mathematical, Physical and Engineering Sciences, 
1897 (1998). 

[2] A. Priel, J. A. Tuszynski, and H. F. Cantiello, in The 
Emerging Physics of Consciousness (Springer, 2006) pp. 
293-325. 

[3] D. Debanne, Nature Reviews Neuroscience, 5, 304 
(2004). 

[4] A. Priel, J. A. Tuszynski, and N. J. Woolf, Journal of 
biological physics, 36, 3 (2010). 

[5] L. Jaeken, IUBMB life, 59, 127 (2007). 

[6] H. F. Cantiello, Journal of Experimental Zoology, 279, 
425 (1997). 

[7] C. Dillon and Y. Goda, Annu. Rev. Neurosci., 28, 25 
(2005). 

[8] L. A. Cingolani and Y. Goda, Nature Reviews Neuro¬ 
science, 9, 344 (2008). 

[9] A. Priel, J. A. Tuszynski, and H. F. Cantiello, Electro¬ 
magnetic Biology and Medicine, 24, 221 (2005). 

[10] A. Adamatzky and R. Mayne, arXiv preprint 
arXiv: 1408.3676 (2014). 

[11] S. Siccardi and A. Adamatzky, Nano Communication 
Networks, 6, 15 (2015). 

[12] S. Siccardi and A. Adamatzky, Submitted (2015). 

[13] A. Adamatzky, International Journal of Modern Physics 
C, 11 , 1321 (2000). 

[14] J. Tuszynski, S. Portet, J. Dixon, C. Luxford, and 


H. Cantiello, Biophysical journal, 86, 1890 (2004). 

[15] E. C. Lin and H. Cantiello, Biophysical journal, 65, 1371 
(1993). 

[16] Z.-s. Ma, J. Wang, and H. Guo, Physical Review B, 59, 
7575 (1999). 

[17] B. Wang, X. Zhao, J. Wang, and H. Guo, Applied 
physics letters, 74, 2887 (1999). 

[18] F. Nakamura, T. M. Osborn, C. A. Hartemink, J. H. 
Hartwig, and T. P. Stossel, The Journal of cell biology, 
179, 1011 (2007). 

[19] H. N. Higgs and T. D. Pollard, Annual review of bio¬ 
chemistry, 70, 649 (2001). 

[20] L. J. Sampson, M. L. Ley land, and C. Dart, Journal of 
Biological Chemistry, 278, 41988 (2003). 

[21] T. Schubert and A. Akopian, Neuroscience, 125, 583 
(2004). 

[22] A. Priel, A. J. Ramos, J. A. Tuszynski, and H. F. 
Cantiello, Biophysical journal, 90, 4639 (2006). 

[23] L. M. Griffith and T. D. Pollard, The Journal of cell 
biology, 78, 958 (1978). 

[24] L. Dehmelt and S. Halpain, Journal of neurobiology, 58, 
18 (2004). 

[25] T. J. Craddock, J. A. Tuszynski, and S. Hameroff, PLoS 
Comput Biol, 8, el002421 (2012). 

[26] C. Joachim, J. Gimzewski, and A. Aviram, Nature, 408, 
541 (2000). 



